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Abstract 



A generalized Strangeness-incorporating Statistical Bootstrap Model (SSBM) 
^ ! is constructed so as to include indepedent fugacities for up and down quarks. 

I Such an extension is crucial for the confrontation of multiparticle data emerg- 

ing from heavy ion collisions, wherein isospin symmetry is not satisfied. Two 
constraints, in addition to the presence of a critical surface which sets the 
boundaries of the hadronic world, enter the extended model. An analysis 
pertaining to produced particle multiplicities and ratios is performed for the 
S + Ag interaction at 200 GeV/nucleon. The resulting evaluation, concern- 
ing the location of the source of the produced system, is slightly in favor the 
source being outside the hadronic domain. 



PACS numbers: 25.75.Dw, 12.40.Ee, 12.38.Mh, 05.70.Ce 



1. Introduction 

Multiparticle production in high energy colhsions is a subject of intense research interest, 
whose history goes almost as far back as that of the strong interaction itself. Indeed, it 
registers as one of the key features entering the analysis of collision processes, involving the 
strong force, at both the experimental and the theoretical fronts. With specific reference to 
relativistic heavy ion collisions, the task of accounting for the produced multiparticle system 
is by far the most important issue to consider for extracting information of physical interest. 

A notably successful theoretical approach, through which experimentally observed parti- 
cle multiplicities have been confronted, is based on the idea of thermalization. Within such 
a context, one views the multiparticle system, emerging from a given high energy collision, 
as being comprised of a large enough number of particles to be describable in terms of a 
thermodynamical set of variables. Relevant, standard treatments appearing in the literature 
[1-5] adopt an "Ideal Hadron Gas" (IHG) scheme, wherein any notion of interaction is totally 
absentQ. 

The fact that such analyses, ranging from e~^e~ to A + A collisions, produce very satis- 
factory results simply verifies, a posteriori, that the thermalization assumption is justifiable. 
Beyond this realization, however, no fundamental insight and/or information is gained with 
respect to the dynamics operating during the process, which produced the multiparticle sys- 
tem in the first place. Given, in particular, that the object of true interest, in the case of 
relativistic heavy ion collisions, is whether the original thermal source of the multiparticle 
system is traceable to a region that belongs, or not, to the hadronic phase, an IHG-type of 
analysis renders itself totally inadequate. Clearly, only if interactions are taken into con- 
sideration does it become relevant to ask whether or not a change of phase has taken place 
during the dynamical development of the system. 

In a recent series of papers [6-8], we have pursued a line of investigation which, on the one 
hand, approaches the study of multiparticle systems from the hadronic side (just as the IHG 



*Only the repulsive form of interaction has been introduced in some cases through a hard sphere 
model. 
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case) while, on the other, incorporates the effects of interactions in a self-consistent way. We 
are referring to the employment of a scheme, known as Statistical Bootstrap Model (SBM), 
which was originally introduced by Hagedorn [9-11], much before QCD was conceived and 
was subsequently developed via notable contributions by a number of authors. Excellent 
reviews articles on the SBM can be found in Refs [12]. 

The crucial feature of the SBM is that it adopts a statistical-thermodynamical mode of 
description, which admits interactions among its relativistic constituent particles via a boot- 
strap logic. According to the SBM, the constitution of the system is viewed at different levels 
of organization (fireballs) with each given level being generated as a result of interactions 
operating at the preceeding one. The remarkable feature of the SBM is that the so-called 
bootstrap equation (BE), which results from the aformentioned reasoning, defines a critical 
surface in the space of the thermodynamical variables, which sets an upper bound to the 
world of hadrons and implies, under precisely specified conditions, the existence of a new 
phase of matter beyond. 

Let us briefly review this construction while giving, at the same time, an overview account 
of the bootstrap scheme itself. We start by displaying the generic bootstrap equation, whose 
final form reads 



where </? is the so-called input function, since its specification involves an input from all 
observable hadrons and G incorporates, via the bootstrap logic, the mass-spectrum of the 

system given in terms of fireballs of increasing complexity. Note that, (/? and G as functions 
depend on a thermodynamic set of variables (temperature and fugacities). 
The key feature of the BE is that it exhibits a square root branch point at 



which defines a critical surface in the space of thermodynamical variables that sets the 
limits of the hadronic phase, in the sense that eq. (1) does not posses physically meaningful 
solutions beyond this surface. This is not to say that the BE is thermodynamically consistent 
with the existence of a different phase on the other side. In this connection, the deciding 



cp{T, {A}) = 2G{T, {A}) - exp(r, {A}) + 1 , 



(1) 



(^(r,{A}) = ln4-l, 



(2) 
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factor is the form of the so-called spectrum function p{m?) entering the definition of G and, 
in particular, the way it factorizes into a kinematical and a dynamical part (see following 
section). 

The final ingredient of the SBM is the employment of a grand canonical partition function 
Z{V,T, {A}), which accounts for thermodynamical properties. In combination with the BE, 
it furnishes a thermal description of a system comprised of relativistic entities (hadrons/ 
fireballs) interacting with each other. 

Recognizing the importance of the role played by the quantum number of strangeness 
in providing possible signals for a presumed QGP phase, we have extended the SBM by 
introducing a fugacity variable for strangeness into the scheme [6,7][|. We shall be referring 
to the resulting extended construction as the "SSBM" . 

Imposing the condition < S* >= we proceeded to study thermodynamical properties of 
the SSBM. Central emphasis was placed on the choice of the spectrum function, in order to 
acheive an acceptable thermodynamical description, consistent with the existence of a phase 
beyond the hadronic one. The end result is encoded into the following relation expressing 
the partition function in terms of the "bootstrap function" G [6,7] 

In ZiV, T, {A}) = ^ f {X})dy , (3) 

where Hq = ,-2^)^34^ with B the MIT bag constant. 

In combination with the critical surface condition furnished by the BE, one is able to relate 
the critical temperature Tq at vanishing chemical potential with B. This occurence provides 
a direct connection between QCD-inspired phenomenology and critical temperature for the 
hadronic state of matter. Our numerical studies have been based on the choice Tq = 183 
MeV, which corresponds to the maximum acceptable value for B^^^, namely 235 MeV. Such 
a choice is consistent with the strangeness chemical potential fig remaining positive definite 
throughout the hadronic phase while maximally extending the region of the hadronic phase 
and thereby rendering our appraisal of the proximity of the source of the multiparticle system 



tThe latter was not included in the original SBM as the main preocupation at the time referred 
to nuclear matter. 
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to the critical surface (or beyond) as conservative as possible. 

Subsequently, we generalized the SSBM [8] by introducing a further "fugacity" variable 
7s, which allows for partial strangeness chemical equilibrium. This extension of the model 
enables us to confront the data with an open perspective on strange particle production, as we 
let the observed particle multiplicities and ratios determine whether strangeness saturation 
has taken place or not. 

In [8] we also conducted a systematic study of multiparticle states (particle multiplicities 
and particle ratios) produced in S* + S* as well as in p + p collisions at CERN (experiments 
NA35 and UA5, respectively), the latter considered more as a test case. Our results yield 
an excellent account of particle multiplicities and ratios (equally good, if not slightly better 
than IHG results [5,14,15]). More importantly, we have identified a region in the space 
of thermodynamical parameters where the source of the produced multiparticle state is 
expected to lie and appraised its location with respect to the limiting surface of the hadronic 
phase. 

For the S + S interaction we summarize the highlights of our findings as follows: 

(a) The quahty of our results were similar to that given by IHG analysis. This further 
justifies the thermalization hypothesis. 

(b) Almost full saturation of strangeness was observed, which accounts for an enhanced 
production of strange particles relative to non-strange ones. 

(c) The source of the multiparticle system was found to he just outside the hmits of the 
hadronic phase, as established by the SSBM. 

(d) An excess of pions (SSBM/experimental = 0.73) is observed, which is not fully com- 
patible with the theoretical prediction of a purely hadronic phase. At the same time, entropy 
considerations also give SSBM/QGP — 0.71 — 0.78, pointing, together with (c), to a source 
being in the doorway of a deconfined phase. 

These findings strongly suggest that in the S+S interaction at 200 GeV/nucleon the ther- 
malized, strangeness-saturated source of the multiparticle system has exceeded the hadronic 
sector and has entered the lower limits of the QGP phase. 

Now, the S + S colliding system is symmetric under isospin transformations, hence con- 
sistent with the simplification = adopted in our previous work on the SSBM. In the 
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present paper we shall further extend the model so as to accomodate isospin non-symmetric 
systems. Such a step will enable us to confront multiparticle data for the S + Ag collision 
experiment (NA35), at CERN. As we shall see, this further extension imposes a new con- 
straint on the system which relates charge and baryon numbers. It follows that the SSBM 
extension we shall be discussing amounts, at the hadronic level, to introducing a fugacity 
variable pertaining to total charge. 

The presentation of the new extension of the SSBM, accomodating isospin non-symmetric 
systems, will be accomplished in Section 2. The profile of the relevant construction, accen- 
tuated by the presence of the critical surface in the space of thermodynamical parameters 
as well as the two surfaces resulting by the imposition of the two physical constraints, will 
be discussed. 

Our confrontation of the data (particle multiplicities and ratios) for the S+Ag experiment 
(NA35 at CERN) is presented in Section 3. Our concluding remarks arc made in Section 
4. Two appendices are devoted to corresponding discussion of a more specialized nature. 
In Appendix A we establish that the value for the critical temperature for zero chemical 
potentials corresponds to a maximum on the critical surface. Appendix B discusses the 
subtle points involved in the minimization of the x^- variable given the presence of constraints 
and the critical surface beyond which the SSBM has no analytical validity. 



2. Isospin non-symmetric SSBM 

In this section we shall realize the construction of a maximally extended SSBM, acco- 
modating both partial strangeness saturation and isospin asymmetry. In this way we shall 
be in position to perform thermal analyses pertaining to (non) strange particle production 
in nucleus-nucleus collisions in which the total number of participating protons differs from 
that of neutrons. 

2a. Preliminaries 

The set of variables in terms of which the initial quantification of the bootstrap scheme is 
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accomplished naturally associates itself with input particle (and fireball) attributes. These 
are number densities and four-momenta pertaining to particle/fireball species. As hinted to 
in the introduction the situation we wish to consider in this paper involves, the following 
number densities: Baryon number b, net strangeness s, overall strangeness \s\ and "net 
charge" q. The employed sequence respects "historical" order in the following sense. In the 
original SBM only b enters, the SSBM construction of Ref. [6,7] includes s(=strangeness 
minus anti-strangeness number) while the extension of Ref [8] has added |s|(=strangeness 
plus anti-strangeness number) to the list. 

Our present effort amounts to a further extension of the SSBM through which we incor- 
porate a "net charge" number density into the bootstrap scheme. To quantify our consider- 
ations regarding this new variable let us focus on the initial states entering a nucleus-nucleus 
collision process and consider the ratio j^, where Np^{N^) denotes the total number of 
protons (neutrons) participating in the collision. Suppose this ratio is equal to unity. It then, 
follows that 

^ ^ 1 (4) 

A^^" < S - < (5 ' ^ ^ 

where < Q and < B are the incoming total charge and baryon numbers, respectively. 
Equivalently, the above condition reads 

<B>^2<Q> . (5) 

By introducing a "net charge" particle density into the bootstrap scheme we declare our 
intention to confront A + A collision processes which do not, a priori, respect the condition 
given by (5). It is not hard to see that the latter corresponds to an isospin non-symmetric 
system, at least as far as its initial (incoming) composition is concerned. 

With reference to (5) the quantification of isospin asymmetry can be parametrized as 
follows 

<B>=(32<Q> , (6) 

where 

ivrin _|_ jvrin 

(3^ P ^ . (7) 

p 
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Our actual preoccupation, of course, is with the description of the produced, final states. 
Accordingly, we shall eventually impose (6) as a constraint on the system. 

We close this general exposition with a brief discussion of the particular version of the 
SSBM we have adopted throughout our work as far as the issue of dynamics vs kinematics is 
concerned. Generically speaking, the SSBM construction involves a mass spectrum function 
T whose dependence is on the set of variables {p^, 6, s, g}. A kinematic factor -B(p^) enters 
the equation (see the BE in the following subsection), the specific choice of which classifies 
different versions of bootstrap models, according to asymptotic behaviour, as the fireball 
mass goes to infinity. Our specific commitment to the form -B(p^) has been discussed at 
great lenght in Refs. [6,7]. We have argued that there are desicive physical advantages in 
favor of the choice 

B(P') = Bip') = 1^ , (8) 

where is the (boosted) four- volume associated with a given particle/fireball and the 
corresponding four-momentum. The two four-vectors being parallel to each other imply a 
relation of the form 

V„ = ^P„, (9) 

V being the rest frame volume. We, therefore, have 

B{p-) B{m^) = ^ . (10) 

The mass mass spectrum aquires the asymptotic form 

f(m^{A})"^C'({A})m-i-"exp[m/r*({A})] . (11) 

The above relations determine the version [a = 4) of the bootstrap model we have found 
to be physically relevant. It should be pointed out that the bulk of the work surrounding 
the bootstrap model, prior to the introduction of strangeness, was based on the choice a — 2 
[13]. 



2b. Construction of the model 

The initial form of the BE reads 
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B{p')fip^b,q,sAs\) = g,,slslHp'W-ml,^^^^ S^p-^^Pi)' 

« n=2 1=1 

input term 

n n n n 

■ ^KiP -Y.^i)Y. ^K{q -Y.ii)12^K{s-Yl Si) Y.^K{b-Yl \s\i) 

bi i=l bi 1=1 Si i=i 1^1^ i=i 

n 

n B{Pi)f{pi, bi, Qi, Si, \s\i)d'^pi . (12) 

i=l 

The new feature, with respect to our previous extensions of the bootstrap model, is the 
introduction of electric charge Q as an additional variablefj. 

Performing five Laplace transforms (one continuous and four discrete) leads to the fol- 
lowing replacement of variables 

(p^ 6, s, \s\,q) {T, Ab, Xs, X\s\,Xq) , (13) 

where the A's represent fugacity variables corresponding to number densities and T is the 
temperature, as recorded in the center of mass frame. 

As the final states are composed of hadrons, rather than just baryons, we find it more 
convenient to pass from the original set of fugacities into one given in terms of valence quark 
fugacities. The transcription is made according to the relations 

Xb = AuA^, Xq = XuX^^, X\s\ = 7s) = ^dXs ^ ■ (14) 

The important implication of the above relations is that they facilitate a thermodynamical 
description of the system in terms of (valence) quark fugacities, thereby enabling us to 
accomodate the presence of any kind of hadronic particle in the final system. Specifically, 
the form of the functions if and G entering the bootstrap scheme is given by 

ip{T,Xu,Xd,Xs,'ys;Ho) = 2'KHQTY,K{Xu,Xd,Xs,'ys)Y.9^imliKi l^j (15) 

a i \ 1 / 

and 



^In our thermodynamical context the electric charge will enter as charge number density q. 
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roo 

G{T, A„, Xd, A„ 7,; Ho) = 2TrHoT / m^To{m\ A„, A^, A„ js)Ki{m/T)dm'' , (16) 

Jo 

where Ki denotes the modified Bessel function of the second kind and where the general 
form of the fugacities Aa, pertaining to the totality of hadronic famihes, is 

Aa({A}) = A^-"-A^'^-"^"A^-"^"7,"^+'^^- , (17) 

where rii is the number of the i quarks contained in the hadron of the a family. For the 
particular case of the fugacities of light unflavored mesons, one can employ the parametriza- 
cion ci{uu + dd) + C2SS, with Ci + C2 = 1, see Ref [8], whereupon the corresponding variables 
assume the form 

Aa({A}) = Ci + C27.'. (18) 

The bootstrap equation (1), written analytically for the case in hand, reads 

ip{T, Xu, Xd, As, 7s) = 2G{T, A„, A^, A^, 7^) - exp[G(r, A„, A^, A^, 7^)] + 1 , (19) 
while the critical surface is determined by (either one of) the relations 

(p{Tcr,Hucr:IJ'dcr:IJ'scrnscr;Ho) = In 4 — 1 (20) 

and 

G{T(~f, Hu cr-i l^d cr-i l^s cri Is cr 'i Hq) = In 2 (21) 

Clearly, the critical surface corresponds to a 4-dimensional surface immersed in the space of 
the 5 thermodynamical variables. 

The constant parameter Hq, related directly to the MIT-bag constant (see remark fol- 
lowing eq (3) and Refs [6,7]), can also be linked to the critical temperature at vanishing 
chemical potentials by 

<^(To, Hu = 0, /Xrf = 0, Us = 0, 7.; i/o) = ln4 - 1 . (22) 

Through this relation Hq can be directly related to To, for a fixed value of the "fugacity" 
7s. At the same time we demonstrate, in Appendix A, that Tq corresponds to the maximum 
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value for the temperature on the critical surface, irrespective of the value 75. In our previous 
work this feature was simply assumed. 

In order to acquire a concrete sense concerning the profile of the critical surface we have 
conducted a number of numerical studies which are displayed in Figs. 1-3. In these figures 
we present various sections of the critical surface, having chosen Hq such that Tq — 183 MeV 
for 7s = 1. 

Fig. 1 depicts projections of the critical surface on the T) -plane for three repre- 
sentative values of fid and three for fig- One observes that the critical surface "shrinks" 
(equivalently, "narrows") as Hs reaches higher positive values, starting from zero. This 
"shrinkage" is more pronounced in the vicinity of vanishing 

Fig. 2 displays critical surface projections on the //(i)-plane for three different values 
of T and Hs- One notices that the projections are (approximately) symmetric with respect 
to the line /i„ = fid- A second point is that the lowering of fig causes an expansion of the 
region occupied by the hadronic phase in the {fiu, fJ'd) plane. 

Finally, Fig. 3 shows projections on the (T, fis) plane of the critical surface for fixed 
values of and of fj,d- We notice that for fixed fid and fj,s the critical temperature falls with 
increasing (absolute) values of The same holds true under the exchange fiu ^ fid- 

2c. Imposition of constraints 

Given the constitution of the initial colliding states, we must impose the constraints 
< 5* >= and < B > —(32 < Q >— Q on the system as a whole. To this end we must refer 
to the partition function for our chosen version of the bootstrap scheme, as given by eq. (3). 
The constraints have the generic form 




(23) 



with 



Fr{y,{\}) = \ 



■s 



(24) 



for the imposition of < 5" >= and 



Hv. {A}) 



1-4/3 , dv{y,{\}) , 1 + 2/3 , d^{y, {A}) 

7^ At, 1 Ad ^TT 



(25) 
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for securing the constraint < B > —(32 < Q >= 0. 

These conditions constitute a system of two equations whose solution yields a 3- 
dimensional hypersurface in the space of thermodynamical variables on which the given 
system is constrained to exist. Let us denote this surface by Hph, where "p/i" stands for 
physical. Clearly, the intersection between Hph and the critical surface defines the limits of 
the hadronic world for the system with the given constraints. This intersection comprises a 
two-dimensional surface whose numerical study is presented in Figs. 6-8. 

Figures 4 and 5 give corresponding perspectives of the profile of Hph whose basic aim is 
to display its variation with p. We have considered the cases /? = 1 (A^^^" = N^), (3 = 2 
(A/";" < iV;") and (3 = 1/2 (A/"™ > A/";")! Fig. 4 depicts projections of Hph in the (T,/i,)- 
plane for fixed values of \u and 7^, while Fig. 5 shows corresponding projections in the 
(a^u, ^^d) plane. From the first figure we record the tendency of ^ig to increase with (3, for 
fixed values of (T, Au,7s). From the second we witness the (expected) behavior /i„ = for 
/? = 1, /i^j < /irf for /5 > 1 and /i^ < /i^ for (3 <1. 

Finally, in Figs. 6-8 we present results of numerical studies pertaining to the intersection 
between Hph and the critical surface. In Fig. 6 the 2-dimensional intersection is projected on 
the T)-plane, for our three representative values of (3. As one might expect, an increase 
of (3 induces a decrease of Hucr for constant temperature. Fig. 7 shows the corresponding 
projections on the /id)-plane exhibiting similar connections between /5-values and the 
relation among ^ucr and ^idcr- In Fig- 8 we consider projections in the (yUu, yUs)-plane. Here 
we surmise that for fixed value of /i„ cr an upward move of (3 with respect to 1 {N^ < N^) 
induces an increase in the (critical) chemical potential of the strange quark. 

This concludes our discussion of the isospin non-symmetric SSBM. We shall proceed, in 
the next section, to confront experimental data encoded in the multiparticle system produced 
in A -|- A collisions, in which we do not have isospin symmetry. 



^For purposes of comparison we have also drawn corresponding projections for the IHG model. 
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3. Analysis of S + Ag data of NA35 

In this section we shall perform a data analysis referring to particle multiplicities recorded 
in the NA35 S + Ag experiment at 200 GeV/nucleon at CERN. The method we shall use is 
similar to the one presented in [8] . The main differences are that our space is described by the 
set of the six thermodynamical variables {VT^/An^, T, {A}), i.e. one more variable is present 
and that the system is subject to two constaints, namely < S >= and < B >= 2/3 < Q >, 
instead of one. The latter will be enforced via the introduction of corresponding Lagrange 
multipliers. 

The theoretical values of the thermodynamical parameters are adjusted via a x^-fit by 
minimizing the function 



x'(VTyAn',T,{X},{l})^Y: 



i=l 



+ ^hHk{T,{X}) . (26) 

k=l 

where 1^ are Lagrange multipliers accompanying the corresponding constraints as given by 
(23) and the A^f ^"^J/ are given by 

^3 /A^3 



jytheory ^ L d\nZ{VTyAn',T,{X},...,K,...) 
V dXi 



(27) 

=Ai=...=l 



The minimization of ammounts to solving the following system of eight equations 

axHx, .,) ^^ (i = l,...,8), (28) 

dxi 

with {x,} = (y^V47r^T,{A},{^}). 

An outline of the procedure involved in realizing a numerical solution of the minimization 
problem has been given in Ref. 10. There, we have also discussed the methodology by which 
we determine correction factors for Bose/Fermi statistics. We shall not repeat the general 
argumentation here, nevertheless we do present in Appendix B a discussion of some technical 
aspects involved in the relevant procedure for the case in hand. 

Turning our attention to the S + Ag collision at an energy of 200 GeV/nucleon, using 
the methodology that has just been described, we set as our first task to specify the value 
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of the /9-parameter appropriate for the process under study. As far as the ^^S nucleus is 
concerned, the single isotope with nucleon number 32 {Z = 16) is employed, whereas for 
silver there are two stable isotopes with nucleon numbers 107 and 109 {Z = 47) entering, 
respectively, a mixture composed of 51.84% and 48.16% fractions. This accounts for an 
average nucleon number of 107.96. It turns out that it makes little difference whether one 
assumes that all the nucleons entering the 5* + Ag system participate in the collision process, 
or that the "active" part of the Ag nucleus is determined by some "realistically" assumed 
geometrical configuration. For our numerical applications we shall fix the value of (3 at 1.10. 
The emerging results are displayed in a series of Tables and Figures. 

In table 1 we present adjusted sets of values for the thermodynamical parameters with 
a corresponding estimation for x^/d^f^. The presented numbers correspond to evaluations 
where all particle multiplicities are taken into account (1st row) and where, in turn, one of 
the particle species is excepted. One notices a decive improvement when (mostly pions) 
are excluded from the fit. This occurence makes meaningful the separate treatment of the 
full multiplicity analysis from the one(s) where pions are excluded. 

The experimental data pertaining to particle multiplicities have been taken from [16- 
20,14,5] and are entered in the first column of Table 2. The second column gives the theoret- 
ical estimates of populations based on the corressponding adjusted set of thermodynamical 
parameters with all particle species included. The third column pertains to the adjusted set 
with the absence of pions. The last column corresponds to the same situation but with the 
critical surface pushed slightly outwards by setting Tg = 183.5 MeV at 7^ = 1. 

Table 3 exhibits the correction factors due to Bose/Fermi statistics for each particle 
species. We have covered each of the three cases entering the previous table: All particle 
species, exclusion of pions with Tq = 183 MeV at 7^ = 1 and Tq = 183.5 MeV at 7^ = 1, 
respectively. Table 4 summarizes the adjustment of the thermodynamical parameters accord- 
ing to the fit (along with the estimate for x^/dof) for each of the three aformentioned 
cases. Finally, Table 5 presents particle ratios (used only for the case where all the multiplic- 



**To ensure that the estimate is carried out without leaving the domain of analyticity of the 
SSBM the numerical values in table 1 correspond to Tq = 190 MeV for 7^ = 1. 
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ities are included) , taken with respect to negative hadron population which has the smallest 
experimental uncertainty (see Ref [7] for a relevant comment). 

Pictorial representation of results with physical significance is given in Figs. 9-12. In 
the first of these figures we display bands, per particle ratio, in the (//.„, T)-plane with 7^ 
fixed at 0.67 (see Table 4). These bands are determined by the experimental uncertainty 
per ratio. The bold solid line marks the boundary of the hadronic world, beyond which the 
SSBM does not present analytic solutions. All particle ratios are used (as per Table 5), i.e. 
pions have not been excluded in the plot. No overlap region of the various bands is observed, 
nevertheless we have marked with a cross the center of a region of "optimum overlap" which 
lies inside the hadronic world. 

Fig. 10 considers corresponding bands of particle populations. Since the variable 
VT^/Att^ also enters our considerations we fix it according to its adjusted value of 1.23 
(sec second column of Table 4). It should be pointed out that the upper limit for the ex- 
perimental Kg population (17) as well as the whole band of the negative hadron population 
(175-197) correspond to fitted values for the thermodynamical parameters that are outside 
the hadronic domain (bold sohd line) . The dotted fine marks the smallest value of on the 
boundary surface. We have determined that 51.6% of particle multiplicities are compatible 
with a source outside the hadronic domain. A zoom around the vicinity of the minimal 
X^- value on the critical surface is presented in Fig. 11. 

A comparison between measured and theoretically determined, according to our fit, 
multiplicities is summarized in Fig. 12. Experimental points, with error bars, are represented 
by heavy dots. Theoretically determined points are marked according to the three cases 
studied throughout, i.e. all particle species inclusion and exclusion of pions with Tq set, 
respectively, at 183 MeV and 183.5 MeV for 7^ — 1. Once again we notice a dramatic 
improvement of the fits when pions {h~) are excluded. As we have argued in Ref. 10, this 
occurence seems to signify an excess of pion production, incompatible with a pure hadronic 
phase (SSBM/experimental= 0.69 ±0.04). 
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4. Conclusions 

In this work we have apphed the SSBM to analyze experimental data from the S + Ag 
collision at 200 GeV per nucleon pertaining to produced particle multiplicities {An projection) 
recorded by the NA35 collaboration at CERN. The quintessential aspect of the model is 
that it accomodates interactions in a self-consistent way within the framework of a thermal 
description of the (relativistic) multiparticle system. Moreover, it designates a precicely 
defined boundary for its applicability, a feature which plays a central role in the assesment 
of our results. It should finally be reminded that for the construction of this bounbary we 
have chosen the largest possible, physically meaningful values of Tq{B) so as to enlarge the 
hadron gas domain and avoid over-optimistic interpretations, regarding the location of the 
source with respect to this boundary. 

Our primary objective has been to locate the source of the multipaticle system in the 
space of the relevant thcrmodynamical set of parameters. In this connection, we have found 
that the data point towards a thermal source that lies just outside the hadronic phase. The 
situation is not as pronounced as in the previously analysed case oi S + S collisions [8] at 
200 GeV per nucleon, where a much larger weight in favor of the source being outside the 
hadronic boundaries was determined. The overall situation resulting from our data analyses 
{p+p, S + S, S + Ag) is depicted in Fig. 13. One notices the proximity of the source location 
for the two nucleus-nucleus collision proccess just beyond the hadron phase as well as the 
(expected) placement of the source for the p+p collision well inside the hadronic domain. 

The experimental data give a substantial excess of pion (entropy) production compared 
to theoretical predictions. This strongly hints that the source of the emerging multiparticle 
system from the S + Ag collision is in the doorway of the QGP phase. As already pointed 
out in the Introduction, an estimate of the entropy associated with the pionic component of 
the produced system for the S + S collision gives [8] a theoretical to experimental ratio which 
is not compatible with hadronic physics and necessitates a location of the source outside the 
hadronic domain. A similar behavior persists in the present case as well. 

A final result of interest was the observance of a tendency towards strangeness saturation. 
This indicates that the source has acheived almost full thermal and chemical equilibrium, as 
expected and required for a phase transition to QGP. 
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We thus conclude that in the S + S and S + Ag interactions at 200 GeV/nucleon we have 
witnessed for the first time the appearance of definite signals linking these interactions with 
the QGP phase. 

The fully extended SSBM is now in position to confront multiparticle data emerging from 
any A + A collision experiment, including strange particles. In this respect, the methodology 
can be applied to other ongoing experiments, e.g. Pb + Pb, as data becomes available and, 
more importantly, on the future ones from RHIC and LHC. On the theoretical side, it would 

be extremely interesting to connect a scheme such as the SSBM coming from the hadronic 
side to corresponding microscopic-oriented accounts of QGP physics [21]. 
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Appendix A 

We shall show that the critical temperature value Tq, as defined in the text, corresponds 
to its maximum value on the critical surface. 
We start by re-expressing eq. (20) in the form 

r = /({A},7.), (29) 

where wc ignore critical value indications on each variable for notational simplicity. 
A maximum for T corresponds to an extremum 

9/({A},7.) 



-0, i=l,...,4. (30) 



dXi 

As long as one remains on the critical surface the above condition can be easily transcribed 



to 



5/({A},7s 



and since Xidip/dT 7^ it must be that 



dif/dXi Xidtp/dXi 
^^~dip/dT ^~ Xidip/dT 



, z-l,...,4 (31) 



'dX, 

Now, for the hadronic fugacities we may write 



^ilr^O' z = l,...,4. (32) 



Aa = ci + A4A}7f- , (33) 

where Af,, ci and A^^a can be read from (17) and (18). (For examble, in the case of the A 
Baryons we have Ci = 0, A^ = A^A^A^ and Ngi, = 1.) Therefore 

dX 

Xi^ = Ar,,A,{A}7f- , ^ = 1, 2, 3 (34) 
oXi 

and 

Xi^ = Ar,,A,{A}7f- , i = 4 . (35) 

Given the above set of equations there will be a corresponding family of antiparticles for 
which 
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dX 

Kg^ = -iV,aA4A}-Si^- , ^ = 1, 2, 3 (36) 

and 

= Ns.h{\}-'l^- , i = 4 . (37) 

will hold true. 

The last four equations applied to (32) give 

E(^'>{^} - h{\Y')Ni,l^'-F,{V, T) = , i = 1, 2, 3 , (38) 

a 

Y.{\{\} + \{\}-')N,,^^-F,{V, T) = , z = 4 , (39) 

a 

where the index "a" runs solely over particles. Next we see that the equation (39) is not 
possible to hold true because the left part is always positive (the numbers A^^a are physical). 
Therefore an extremum of the temperature with respect to 7^ does not exist. Turning to 
(38) we observe that a solution could be found if for all {A} we had 

A6{A} - \,{\}-^ = 0, V6 4^ \{\} = 1, V6 . (40) 

An obvious solution for (40) is 

A« = Ad = As = 1 <^ /i„ = /^d = /^s = 1 ■ (41) 

The last equation defines an extremum for the critical temperature with constant 7^. 
On the other hand we have 

^ ^^^Ir) = Y- El-^f {-^i + \W)Nr.N^,^^-F,{V, T) > , z = 1, 2, 3 . (42) 



dXj \ dXiJ Xj a 

That is, for every value of i = 1, 2, 3 each one of the above equations, once two values among 
the {A} are fixed, will have a unique solution. That happens, because from (42), one can infer 
that Aj^ is a genuine rising function with respect to Xj and so it has a unique solution. By 
extension the simultaneous solution of the three equations will be unique. So the extremum 
we have calculated is unique. This extremum cannot correspond to minimum, since the 
critical surface has zero critical temperature for non zero chemical potentials (e.g. see Figs. 
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1,3). Since always T > 0, if the point which corresponds to (41) was a local minimum, then 
we should have another extremum somewhere else, which is imposible, since the extremum 
is unique. Therefore (41) corresponds to a total maximum for a given value of 7^. 



Appendix B 

In our study we have to calculate r constraints {r — 1 for isospin symmetry and r — 2 
for isospin non-symmetry) and different particle multiplicities as functions of the thermody- 
namical variables (T, {A}). In general all these quatities can be written as 

where 

Rj = Hj,j<r R^ = Nf_r\3>r, (44) 

where H and N^'^'^'^'^ are given from (23) and (27), respectively, and 

VT^ \dip{y,{X},---,Xj~r,---y 



Qj — 1 j ^ f 1 Qj 



, J > r , (45) 
=1 



with Fj given from (24) and (25). 

In order to evaluate the optimized set of variables (T, {A}) in our large working space we 
have to turn to the use of the generalized Newton- Raphson method which converges quickly 
but requires the knowledge of the derivatives of (43). If we try to calculate these derivatives 
with respect to a fugacity Aj from (43) we find 

dR,{T,{\}) 



T 



(T dy j exp[G{y, {A})] d^{y, {A}) 
Jo y^\{2-eMG{y,{m' dX, 



+ 



1 dQ,{y,{X}) \ 
2 - exp[G{y, {A})] ^A, j 



(46) 



With the use of the above equation the Newton-Raphson method can proceed for all 
points of the hadronic space which are not close to the critical surface. Problems, however, 
are encountered when the quantities (46) have to be evaluated near and even more on the 
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critiacal surface. When T — > Tcr the function to be integrated in (46) contains a non- 
integrable singularity of the form (2 — exp[G])~^. This singularity cannot be integrated even 
if we use the variable z = 2 — exp[G'(i/, {A})], as we did in Rcfs [7,8]. 

On the other hand the quantities Rj can be expressed as functions of a new set of variables. 
This new set can be formed if we replace the temperature T in favour of the function (/?, or 
equivalently the z variable. Then the Rj can be given from 



Ji z 



dz 



J z 



2 - exp[G{y, {A})] 



(47) 



y dy 

To be able to proceed with the Newton-Raphson method, in this case, we calculate instead 
of (46), the derivatives of Rj when z is constant, i.e. derivatives of the form ^'^•'g^'.^^^^ ^■ 
These derivatives should not present any singularity for any value of z (even for 2; = when 
we are on the critical surface) because, as can be seen from (47), the Rj can be evaluated 
for any values of z and {A}. 

In order to proceed with the evaluation of the derivatives of (47) we can assume that 
we are standing on a surface of constant (p or equivalently constant z. We then let the 
variation to the fugacity Aj be dXi without ever leaving the above mentioned surface. Then 
the variation in z is 



Since using the BE we have 



dz = ^dXi 
dXi 



dG 

d(f 







we arrive at 



dz 
dXi 



dz dG dip 
dG dip dXi 



3^ dX, 



From the last two equations we conclude that 

dip{z, {A}) 



9A,; 







(48) 



(49) 



(50) 



(51) 



Let us comment that on the critical surface 2 — e'^ = 0, so again equation (51) is to hold if 
equations (48) and (50) are to be fulfilled. If we then express the z variable as function of 

the temperature y and the fugacities, z = z{y, {A}), then equation (51) leads to 
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9<p{y: {A}) _ dy_ 
dy dXi 



I My, {A}) ^ Q ^ dy_ 
dXi dXi 



d(f/dXi 
dip/dy 



(52) 



The last relation shows us how temperature is varied with the fugacity Aj on a surface of 

constant z. 

Using the definition 



vAy,W) 



My, {A}) ' 

^ By 



the derivatives we seek can be expressed as 



dR,{z.{X\) 



dX, 



/■^ dz 

Ji ¥^ 



dX, 



But 



vM^,{M),m _ vAy,W] y , vAy,{>^}] 



dX, 



y K 



+ 



dX, 



(53) 



(54) 



d(p \ dQj 5 f dip 



y ^ 



dy J dy y \ dy 



y \dy J dy 



d(f/dXi 
dip/y 



+ 



+ \y 



,d(p\ ^ dQj Qj ( dip\ ^ d'^(p 



dy J dXi y^ \dy J dydX. 



-5 / '^^ 



-2 r 



Qy , 



dp dQj dip dQj 
dy dXi dXi dy 



+ Q, 



d'^p ^ d'^p dp/dXi dp 



dydXi dy'^ dp/dy ydXi 



Prom the last two equations we conclude that 



dRj{z,{X]) 



dX, 



dz \ _^ I dp^ 
1 tildez — 2 I ^ V ^2/ y 



dp dQj dp dQj 



dy dXi dXi dy 



+ 



(55) 



Qj 



d'^p d'^pdp/dXi 5 dp 



z^2- exp[G(|/, {A})] 



(56) 



'•^ \ dydXi dy"^ dp/dy y dXi 

As it is known the function to be 2;-intergrated has not so good behaviour near z — 1. 
So it is better to break the above integral in two parts: 

^^.(^,{A}) 
dX, 
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^1 dy _5 / dip' 
[dy 



dip dQj 

dy dXi dXi dy 



dip dQj j^Q f ^"^V d'^ipdip/dXi dip 
^ [dydXi dy^ dcp/dy y dXi 



+ 



Jzi Z 



dz 



'dip' 

9y, 



-2 r 



dip dQj dip dQj 



dy dXi dXi dy 



f d'^ip d'^ipdip/dXi 5 dip'' 



2 - exp[G{y, {A})] 



• (57) 



^ _1_ j a,[.(.,{A}),{A}] j- ^^^^^^^ ^^^^^ ^^^^ 



(58) 



'■^ \dydXi dy^ dip/dy ' |/c?Ai 

In the above relation zi — 2 — exp[G'(Ti, {A}) and a good choice is zi — 0.5. If 2; > 0.5 we 
are not close to the critical surface and the second integral does not have to be calculated. 

So in general we can set Zi = max{z, 0.5}. 

The derivatives of Rj with respect to z can be calculated easily. They simply read 

^^.(^,{A}) 

A, 

With the above relations the minimization of the function can proceed with the use 
of the Newton-Raphson method. Relation (58) can also be used to find out whether the 
absolute minimum of is outside or inside the critical surface. Suppose we locate the 
minimum value of = (x^)i on the critical surface and this value corresponds to the point 
{z, {A}) = (0, {A°}). Then the absolute minimum of x^ is located inside the hadronic phase 
if for this point we have 

dx'(z, {X'}) 



dz 



< . 



z^O 



(59) 



If the above relation is not fulfilled the absolute minimum of x^ li^s on the outside. 

An alternative method to verify the same thing consists of locating the minimum value of 
X^ on a surface near the critical one inside the hadronic phase. Let this value be (x^)2- The 
absolute minimum of x^ is located inside the hadronic phase if (x^)2 < (x^)i and outside 
otherwise. For the two fits we have performed for S + Ag we had to process 192 points. All 
these points have given us the same results with the use of the two methods. 
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Table Captions 

Table 1 The fitted parameters and the 'x^/dof values for different fits through SSBM in 
the experimentally measured full phase space multiplicities in the S + Ag interaction. In the 
first fit all the multiplicities are included while in the following fits we exclude each time one 
mutipUcity. Tq is set to 190 MeV. 

Table 2 Experimentally measured full phase space multiplicities in the S + Ag interaction 
and their theoretically fitted values by the SSBM, with the inclusion of the hT multiphcity 
and without it (cases A: Tq = 182.94 MeV and B: Tq = 183.5 MeV). 

Table 3 Calculation of the correction factor fi = W«G-^F)i-W^G-Bo)i f particle 

(NiHG-BO)i ^ 

species measured in An phase space in the S + Ag interaction. For the calculation of fi 
the IHG formalism has been used, while the thermodynamical variables have been extracted 
from the SSBM fit with h~ and without h~ (cases A and B). 

Table 4 Results of the analysis by SSBM of the experimental data from the S + Ag inter- 
action (Att phase space), with the inclusion of the h~ multiplicity and without it (cases A 
and B). 

Table 5 Particle ratios from the experimentally measured full phase space multiplicities for 

the S + Ag interaction used in the analysis with . 
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S + Ag (NA35) Full phase space 



Excluded 


T (MeV) 


K 7s 




xVdof 


none 


170.943 


1.540 1.582 1.088 0.662 


2.883 


10.35/3 




162.411 


1.551 1.588 1.132 0.749 


3.869 


4.72/2 


A 


170.363 


1.523 1.563 1.093 0.616 


3.063 


9.41/2 


A 


158.361 


1.625 1.665 1.174 0.636 


4.239 


7.58/2 


P 


175.180 


1.482 1.520 1.067 0.596 


2.807 


4.55/2 


p-p 


170.626 


1.538 1.579 1.090 0.665 


2.904 


10.33/2 


B-B 


172.040 


1.554 1.597 1.083 0.656 


2.786 


8.52/2 


h- 


180.779 


1.642 1.702 1.011 0.839 


1.261 


1.64/2 



Table 1. 
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S + Ag (NA35) Full phase space 



Particles 


Experimental Calculated Calculated Calculated 
Data with h~ without without /i^ 

(Case A) (Case B) 




15.5 ± 1.5 


17.613 


15.181 


15.155 


A 


15.2 ± 1.2 


14.490 


15.424 


15.429 


A 


2.6 ±0.3 


2.3998 


2.5502 


2.5538 


P 


2.0 ±0.8 


3.4614 


2.3612 


2.3547 


p-p 


43 ±3 


42.600 


40.931 


40.937 


B-B 


90± 10 


101.39 


99.307 


99.325 


h- 


186 ± 11 


170.68^ 


128.85'' 


128.57^ 



A correction factor 1.0236 has been included for the effect of Bose statistics. 

A correction factor 1.0190 has been included for the effect of Bose statistics. This multiplicity is not 
included in the fit. 

A correction factor 1.0188 has been included for the effect of Bose statistics. This multiplicity is not 
included in the fit. 

Table 2. 



S ± Ag (NA35) Full phase space 



Particles (i) 


fi{%) for the 
fit with h" 


fi{%) for the fi{%) for the 
fit without h~ (A) fit without (B) 


Ks' 


0.462 


0.641 


0.645 


A 


0.053 


0.232 


0.237 


A 


-0.142 


-0.410 


-0.418 


P 


0.146 


0.105 


0.105 


p — p 


-0.388 


-0.431 


-0.432 


B-B 


-0.434 


-0.417 


-0.417 


h- 


2.357 


1.897 


1.877 



Table 3. 
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S + Ag (NA35) Full phase space 



Fitted 


Fitted with h 


Fitted without h 


Fitted without h 


Parameters 




(Case A) 


(Case B) 


T (Me\f\ 

-L I IV J.C V j 




176 3 




x 


1 544 ± 046 


1.640 


1 641 ± 074 




1 586 ± 050 


1.700 


1 701 ± 084 




1 084 ± 036 


1.012 


1 Oil ± 047 


IS 


670 ± 073 


0.836 


0.84 ± 0.12 




2.74 ± 0.71 


1.23 


1.22 ± 0.60 


/do f 


9.37 / 3 


1.654 / 2 ^ 


1 652 / 2 


l^u (MeV) 


74.1 ±5.7 


87.3 


87.6 ±8.0 


fid (MeV) 


78.7 ±6.0 


93.5 


93.9 ±8.8 


l^s (MeV) 


13.8 ±5.7 


2.2 


2.0 ±8.1 


PiNSIDE 


100% (128/128) 


48.44% (31/64) 





* It is the minimum of x'^ within the Hadron Gas with Tq = 183 MeV (for 7^ = 1), not the absolute 
minimum. 

Table 4. 



S + Ag (NA35) Full phase space 



Particle ratios for 


Experimental 


the fit with h" 


Values 


KVh- 


0.0833 ±0.0095 


A/h- 


0.0817 ±0.0081 


A/h- 


0.0140 ±0.0018 


p/h- 


0.0108 ± 0.0043 


p — p/h^ 


0.231 ±0.021 


B - B/h- 


0.484 ± 0.061 



Table 5. 



Figure Captions 

Figure 1 Projections on the (/i^, T')-plane of intersections of constant /id and constant /ig 
of the critical surface (p{T, /z^, /i^, 7^) = ln4 — 1 for Tq — 183 MeV at 7s = 1. 

Figure 2 Projections on the (/i^, /i„)-plane of intersections of constant T and constant Hs 
of the critical surface (p{T, /i^, l^'di l^s, Is) = ln4 — 1 for Tq = 183 MeV at 7^ = 1. 

Figure 3 Projections on the (T, //s)-plane of intersections of constant //„ and constant 
of the critical surface (p{T, A*s, 7s) = ln4 — 1 for Tq — 183 MeV at 7^ = 1. 

Figure 4 Projections on the (T, //g) -plane of intersections, at fixed {Hu/T — 0.4), of the 

< 5" >= and < S >= 2/3 < Q > surfaces for the SSBM and the IHG for different values 
of 13. For the SSBM case To is set at 183 MeV. 

Figure 5 Projections on the /irf)-plane of intersections, at fixed A„ {f^u/T = 0.4), of the 

< S >= and < B >= 2(3 < Q > surfaces for the SSBM and the IHG, for different values 
of (3. For the SSBM case Tq is set at 183 MeV. 

Figure 6 Projections on the (T, /i„)-plane of the intersection of the < S >= and < B >= 
2(3 < Q > surfaces for the SSBM with the critical surface for different values of (3. Tq is set 
at 183 MeV. 

Figure 7 Projections on the /irf)-plane of the intersection of the < S >= and < B >= 
2(3 < Q > surfaces for the SSBM with the critical surface for different values of (3. Tq is set 
at 183 MeV. 

Figure 8 Projections on the (yU^, yUs)-plane of the intersection of the < S >= and < B >= 
2(3 < Q > surfaces for the SSBM with the critical surface for different values of (3. Tq is set 
at 183 MeV. 
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Figure 9 Experimental particle ratios in the T)-plane for the S + Ag interaction mea- 
sured in 47r phase space with 7^ set to 0.67. The point and the cross correspond to the 
fit with the . The thick solid line represent the limits of the hadronic phase (HG) as set 
by the SSBM. 

Figure 10 Experimental particle multiplicities in the T)-plane for S + Ag interaction 

measured in Att phase space with 7^ set to 0.84 and VT^/Air^ set to 1.23. The point repre- 
sented by the solid circle corresponds to the location of the least value within the hadron 
gas of x^, without the h~. The hues which correspond to — 17, — 175 and h~ — 197 
lie outside the hadronic domain, as set by the S-SBM. 

Figure 11 The same diagram as Fig. 10, but with an enlargement of a smaller area to 
show the common overlapping region (shaded area) within the hadronic phase which is 
compatible with all the measured multiplicities, except h". The lines wich corespond to 
A = 16.4, A = 2.9, p = 1.2, p — 2.8 and p — p = 46 are outside the region of the diagram 
and enclose the shaded region. 

Figure 12 Comparison between the experimentally measured multiplicities in Air phase 
space and the theoretically calculated values in the fit with h~ and without h~ (cases A and 
B) for the S+Ag interaction. The difference is measured in units of the relevant experimental 
error. 

Figure 13 (/is, T)-phasc diagrammc with points obtained from fits to p [8], S + S [8] 
and S + Ag data and corresponding critical curves given by SSBM. 
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S + Ag , full phase space, =0.84, VT^ /4n^ =123 
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